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Abstract 

Structural and electronic properties of hypothetical zinc blende TPGai_a;N 
alloys have been investigated from hrst principles. The structural relaxation, 
preformed within the LDA approach, leads to a linear dependence of the 
lattice parameter a on the T1 content x. In turn, band structures obtained 
by MBJLDA calculations are significantly different from the corresponding 
LDA results. The decrease of the band-gap in Tl-doped GaN materials (for 
X < 0.25) is predicted to be a linear function of x, i.e. 0.08 eV per atomic % of 
thallium. The semimetallic character is expected for materials with x > 0.5. 
The obtained spin-orbit coupling driven splitting between the heavy-hole and 
split-off band at the T point of the Brillouin zone in Tla;Gai_a,N systems is 
signihcantly weaker when compared to that of Tl-doped InN materials. 

Keywords: nitride semiconductors, doping with thallium, electronic 
structure, band-gap 


1. Introduction 

Semiconducting group-III nitrides are well known materials used for fabri¬ 
cation of optoelectronic devices jlj . Parent compounds adopting the wurtzite 
phase, with direct band-gaps starting from 0.65-0.69 eV for InN 0,3 8, 
through 3.50-3.51 eV for GaN 0,8) to 6.1 eV for AIN 8; allow a design 
of ternary alloys within a wide interval of available wave lengths. A further 
band-gap engineering, leading to the infrared devices, could be realized by a 
doping with elements exhibiting a high atomic number. 

Although some Tl-doped III-V group semiconductors have alrea^ been 
obtained, e.g. TGIni-^^P 0], TbGai_a;As, and Tblni-^jAs alloys 0, 8), the 
synthesis of such nitrides has not been reported yet. Meantime, the structural 
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and electronic properties of hypothetical Tl^-Ali.^-N [ld,lll[, ThGai.^^N [12 


and Tl 2 .Ini_ 2 .N [13|| systems were investigated only from hrst priciples. 


Standard DFT calculations within the local density (LDA) and general¬ 
ized gradient (GGA) approximations for TIN predicted a semimetallic char¬ 
acter of this compound with the zero band-gap BE HQ. The recent 
fully relativistic study within the modihed Becke-Johnson (MBJLDA) ap¬ 
proach suggested that the strong spin-orbit coupling (SOG) in zinc blende 
(ZB) TIN may cause an opening of a narrow gap of 0.11 eV at the T point of 
the Brillouin Zone (BZ) Q], similarly to the case of HgS Q]. Such a feature 
of the band structure is exhibited by topological insulators B, E- 

Despite the fact that the ThGai_a;N alloys have been already investi¬ 
gated with LDA/GGA methods Q] some properties of those systems are 
still unknown. Since the standard DFT-based calculations are insufficient 
for a prediction of band structures for narrow band-gap systems, the values 
of the band-gap as well as the SOG-driven separation between the heavy-hole 
and the split-off band at the F point of the BZ {Aso) require studies with 
more advanced methods. 

In this work, the band structures of ZB Tl 3 .Gai_ 3 ;N alloys are investigated 
with the MBJLDA potential [Q. This well-known method has already been 
examined for semiconductors exhibiting narrow band-gaps, e.g. InAs and 


InSb |2l(], and also used in the former study of TGIni_j.N alloys [13|. The 
main aim of this work is a prediction of the band-gaps of Tl-doped GaN 
materials. 


2. Computational details 

The electronic structure calculations for ZB TLGai .t,N systems have 
been performed with the use of the ABINIT code 22|, [2^. First, the equi¬ 


librium geometries were found via stresses/forces relaxation for PAW atomic 
data sets generated with the Atompaw package B] with the Perdew-Wang 
25| parametrization for the exchange-correlation energy. Then, the calcula¬ 


tions of the band structures with the use of the norm-conserving fully rel- 
ativistic pseudopotentials were performed with the MBJLDA (TB09) func¬ 
tional j20j. The pseudopotentials were generated using the Atomic Pseu¬ 
dopotential Engine (Q. Such a complex calculation scheme was necessary 
to obtain reasonable results of LDA-predicted structural parameters and 
MBJLDA band-gaps with SOG effects included. A supercell consisting of 
16-atoms (2x2x2 multiplicity of primitive fee cell) was used to simulate 
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the alloy. The most possibly uniform arrangements of atoms were chosen. 
The total energy convergence in the plane wave basis was sufficient for 50 
Ha energy cut-off and the 4x4x4 k-point mesh with standard shifts for fee 
lattice. 


3. Results and discussion 


The computed lattice parameter a, band-gap (Eg), and split-off energy 
(A 50 ) for parent compounds, zinc blende GaN and TIN, are gathered in Ta¬ 
ble [U As seen, the structural relaxation with the LDA potential for GaN 
leads to very good agreement between the calculated result and the accesible 
experimental data. In the case of TIN, only theoretical predictions are avail¬ 
able in the literature. The calculated here lattice parameter a of this system 
is close to other LDA results, being smaller than the GGA one. 

It is a well-known fact, that band-gaps of insulators are underestimated 
by the LDA approach. The obtained (LDA) value of the band-gap for GaN, 
Eg = 1.69 eV, is almost two times lower than the experimetal one (3.30 eV 
( 3 ^). The MBJLDA-calculated Eg of 3.04 eV is also slightly too low. It is 
worth noting that this value is very close to the result of the full-potential 
calculations. Eg = 3.03 eV, reported in Ref. 


27[. Therefore, the accuracy of 


the pseudopotential scheme used here seems to be satisfactory. Furthermore, 
the obtained in both approaches (LDA/MBJLDA) value of Aso = 16 meV 
for GaN is in very good agreement with the experimental Aso = 17 meV 
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The MBJLDA band structures for GaN and TIN are presented in Fig. [T] 
One can notice that the overall shape of valence bands is similar in both 
systems. However, the Tl-based compound exhibits a semimetallic character 
with an inverted band-gap 1^, beeing related to the anomalous order of 
N 2p bands, since the quartet N 2p^/2 is below the doublet N 2pi/2 fl3j |. 
The strong SOG in TIN leads to an opening of a narrow gap at the T point 
of BZ, estimated to be equal 0.03 and 0.10 eV from LDA and MBJLDA 
calculations, respectively. The high atomic number of T1 is also reflected in 
significant splittings of heavy-hole bands, e.g. in the X — T line of BZ, as 
can be seen in Fig. [TJ 

Because high values of Aso for TIN, ranging from 1.4 to 2.0 eV, have 


been reported in the literature [13|, \l% llTi, one may consider that the po¬ 


sition of the light-hole bands of such a narrow band-gap system is affected 
by the overestimation of the negative band-gap in the standard LDA/GGA 
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calculations. Thus, the lower MBJLDA value of A 50 ~ 1-5 eV, seems to be 
more adequate. It is worth noting that the obtained Aso is slightly higher 
than the results of former MBJLDA study for the same volume of the unit 
cell [ 1 ^. This discrepancy is related to the fact that in the present study the 
calculations were performed fully ab initio. Namely, the parameter c in the 
MBJLDA exchange-correlation functional was not tuned empirically. 

The calculated lattice parameters a for ZB Tl,j,Gai_ 2 ;N alloys, presented 
in Fig. [2], exhibit a linear behavior (the Vegard’s law), similarly to the 
former results for TGAli_ 3 ;N 1^ and Tla;Ini_ 3 ;N [l^ materials. The strong 
bowing of lattice parameters a as a function of a T1 content x reported in 


the former study for Tla,Gai_a,N [ 12 (] systems may be connected with the 


lack of a relaxation of atomic positions during the volume optimisation of 
the unit cell. Such calculations lead to similar results to those of the virtual 
crystal approximation, which have been discussed for some nitride systems, 
e.g. Ina,Gai_j;N and AlNi_ 3 ;Pj; ^|. 


The dependence of calculated within the MBJLDA approach band-gaps as 
a function of a T1 content in Tla,Gai_ 2 ;N alloys, depicted in Fig. [31 exhibits an 
almost linear behavior for x < 0.25 (for in the range from about 3 to 1 eV). 
An an alog ous linear dependence of Eg in Tl,i.Ini_ 3 ,N systems was reported ab 


initio 13 


although some group-III nitride semiconductors exhibit a strong 
band-gap bowing, e.g. Ina,Gai_j;N and Ina,Ali_ 3 ;N alloys [^. Interestingly, 
the experimental studies for Sc-doped GaN and AIN materials also reported a 
linear behaviour of Eg as a function of a Sc content j^, One can consider 
that the possible solubility of thallium in GaN is rather questionable for high 
contents of this element in such a material. Therefore, the linear reduction of 
Eg Ri 0.08 eV per T1 atomic % was predicted arbitrarily, basing only on data 
for X < 0.25. Since the band-gap of Tlo.sGao.sN is almost zero, the T^.Gai_ 3 ;N 
systems are predicted to be semimetallic for x > 0.5. It is worth noting that 
the presented here dependence of Eg for Tl 2 ;Gai_ 3 ;N is significantly different 
from that reported earlier [l^. Namely, the drop of Eg is less sharp in the 
LDA results, which cannot be corrected by a simple shift of bands (the so- 
called ’scissor operator’). Therefore, the band structure of semiconductors 
doped with heavy elements should be investigated with methods beyond 
the standard (LDA/GGA) DFT approach, e.g. hybrid exchange-correlation 
functionals or the used here MBJLDA potential. 

A strong enhancement of SOG-driven splitting between the heavy-hole 
and split-off band at the F point of BZ for Tl 3 ,Gai_ 3 ;N alloys, related to the 
high atomic number of thallium, is expected to be pronounced. However, 
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the results of fully relativistic MBJLDA calculations indicate a significantly 
weaker influence of T1 atoms on the split-off band in the GaN host systems, 
i.e. Aso < 0.5 eV for x = 0.25, when compared to the case of ThIni_ 3 ;N (l3j |. 
Furthermore, for x > 0.25, the split-off band is even shifted to the vicinity 
of the valence band maximum. This effect may be explained by the fact 
that the semimetallic systems with a higher T1 content exhibit an inverted 
band-gap. Therefore, the arrangement of particular bands in materials being 
close to the inversion of the band-gap may be complex and investigations of 
such an issue are beyond the scope of this work. 

4. Conclusions 

The structural and electronic properties of ZB TGCai-^jN alloys were in¬ 
vestigated within LDA/MBJLDA calculations. The linear dependence of the 
lattice parameter a as a function of a T1 content was revealed. The singnif- 
icant decrease of the band-gap of Tl-doped GaN materials is also predicted 
as a linear function of a T1 content x for x < 0.25, whereas the semimetallic 
character is expected for x > 0.5. The obtained SOG-driven splittings be¬ 
tween the heavy-hole and split-off band at the T point of BZ in Tla;Gai_a;N 
alloys are relatively low when compared to those of Tl-doped InN systems. 
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Table 1: Lattice parameter a, band-gap (Eg), and split-ofF energy (Aso) of zinc blende 


GaN and TIN. 

a (A) Eg (eV) Aso (eV) 


GaN: 

this work LDA 
this work MBJLDA 
LDA Ref. [27] 
MBJLDA Ref. [27] 
LDA Ref. [28] 
exp. Ref. [29] 
exp. Ref. [30] 

TIN: 

this work LDA 
this work MBJLDA 
MBJLDA Ref. [1^ 
LDA Ref. \J£ 
GGARef. [17] 


GaN 



4.481 

1.69 

0.016 

- 

3.04 

0.016 

4.49 

- 

- 

- 

3.03 

0.019 

4.49 

- 

- 

- 

3.30 

0.017 

5.136 

0.03 

2.078 

- 

0.10 

1.499 

- 

0.11 

1.397 

5.133 

- 

2.000 

5.153 

_ 

1.830 


TIN 



Figure 1: Band structures of zinc blende GaN (left panel) and TIN (right panel) calculated 
within MBJLDA. 
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Figure 2: Calculated (LDA) lattice parameters a of zinc blende Tla;Gai_a:N alloys as a 
function of thallium content x. 



Figure 3: Calculated by MBJLDA band-gaps (Eg) and split-off energies (Ago) for 
Tlj^Gai-a^N alloys as a function of thallium content x. 
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